Quantum simulation of artificial Abelian gauge field using nitrogen-vacancy center 
ensembles coupled to superconducting resonators 
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We propose a potentially practical scheme to simulate artificial Abelian gauge field for polaritons 
using a hybrid quantum system consisting of nitrogen- vacancy center ensembles (NVEs) and super- 
conducting transmission line resonators (TLR). In our case, the collective excitations of NVEs play 
the role of bosonic particles, and our multiport device tends to circulate polaritons in a behavior like 
a charged particle in an external magnetic field. We discuss the possibility of identifying signatures 
of the Hofstadter "butterfly" in the optical spectra of the resonators, and analyze the ground state 
crossover for different gauge fields. Our work opens new perspectives in quantum simulation of con- 
densed matter and many-body physics using hybrid spin-ensemble circuit quantum electrodynamics 
system. The experimental feasibility and challenge are justified using currently available technology. 
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Gauge field theory 1] , which was developed for describ- 
ing the subatomic interaction, has deepened our under- 
standing of a wide range of physical phenomena: quan- 
tum magnetoresistance oscillations, superconducting vor- 
tices, quantum Hall effects (QHE), and Hofstadter but- 
terfly (HB) [2J etc. However, these phenomena associ- 
ated with the gauge field are very difficult to be directly 
observed in ordinary solid-state systems due to the re- 
quirement of extremely high field. Recently, following 
the idea of quantum simulation [3], several theoretical 
schemes have been proposed to generate artificial gauge 
field in atomic, molecular, and optical systems. The first 
route is taken in cold atomic systems. The gauge field 
can be generated by rotating the trap or lattice UHl, by 
introducing appropriate phase factors for hopping ampli- 
tudes through combining Raman-laser-assisted tunneling 
and lattice acceleration methods [6l-[l0|. or by tailoring 
orthogonal laser-beam coupled to the atomic degenerate 
internal states [ill, [13] , or by simply using the ordinary 
tunneling in an optical lattice jl3| . The second route 
is based on cavity quantum electrodynamics (QED) sys- 
tems, where the gauge field can also appear through ma- 
nipulating the phase factors of the hopping term, in cir- 
cuit QED cavity arrays [l3 | , in confined ions in microtrap 
array [l3|, in coupled resonator optical waveguides (l6j |. 
and in solid-state photonic structures [l3l- In addition, 
the fractional quantum Hall effect (FQHE) can also be 
simulated in optical-cavity arrays by trapping three-level 
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FIG. 1: (Color online) Schematic circuit for the resonator 
array, where the NVEs are magnetically coupled to the quan- 
tized field of the TLR, and the circles denote the central cou- 
pler. Microwave photons can tunnel between adjacent res- 
onators via capacitive coupling. The lattice sites are repre- 
sented by a{px + qy) with a the spacing of the lattices, where 
the dotted grid denote the 'vertical' and 'horizontal' direction. 



atoms and using elaborated laser driving [18]. However, 
due to the required rigorous conditions, up to now only 
few pioneering experiments have successfully mimicked 
the effect of a light-induced artificial magnetic field in 
an optical lattice |20| and a synthetic electric field Bose- 
Einstein condensates [2l| by using a Raman-assisted tun- 
neling method [§] , and an effective time-dependent vector 
potential, respectively. 

In this paper we propose an alternative theoretical ap- 
proach for quantum simulating an effective gauge po- 
tential in a hybrid solid-state system. We consider 
NVEs [22[j28t which are confined in two-dimensional 
(2£)) square lattices of superconducting microwave TLR 
[29| . The NVEs are driven by series of microwave sources 
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FIG. 2: (Color online) Level structure of a NVE, where the 
electronic ground state is an electron spin triplet state (S=l), 
and Dgs/2ir = 2.87 GHz is the zero-field splitting between the 
rus = sublevel and the = ±1 sublevels in the absence 
of external magnetic field. The degeneracy between states 
ma = ±1 is lifted by applying a static magnetic field with 
40 mT parallel to the chip and along the [111] crystalline 
axis within a few degrees. With this orientation, the four 
equivalent NV center crystalline orientations (111) all make 
approximately the same angle 55° with the magnetic field 
so that their resonance frequencies are approximately equal. 
The transition between \g} and |e) with Larmor frequency lj^ 
is coupled to the two-mode TLR with the coupling rates Qa 
(gt) and detunings Aa (At), and is also driven by a classical 
field with the Rabi frequency Q, and the related phase 9. 



with different phases at different sites. We carefully tune 
the frequencies of the cavity modes, the driving strength 
and the phase of microwave sources, in order to induce 
the non-zero tunneling phases around a closed loop in real 
space and to generate a non-vanishing artificial gauge 
field experienced by the orbital motion of the polari- 
tons in NVEs. Our approach utilizes a space-dependent 
coupling between internal ground states of the NVEs, 
which yields crucial phases to create considerable syn- 
thetic gauge fields. This intriguing behavior is indeed 
the analogue of the motion of charged particle in an real 
magnetic field. 

The main merits of the present system include the 
in situ tunability of the parameters of the circuit el- 
ements, individual addressing, the peculiar characteris- 
tics of the NVE [sO] (e.g., long coherence time at room- 
temperature ), a nd the scalability of cavity resonator ar- 
rays IsTI. l32j. Very recently, D. Underwood et al ex- 
perimentally demonstrated 25 arrays of twelve capaci- 
tively coupled TLRs and accessed the feasibility of quan- 
tum simulation in cavity QED systems 's^. It implies 
that the polariton-polariton interaction between distant 
NVEs can be effectively tuned in a controllable way, 
which renders our scheme to be more practical. By 
combining the related spectroscopic circuit QED tech- 
nique for readout of the quantum states of individual 
constituent elements, we could probe the properties of 
the system by independently detecting the correlation 
between distant sites or output fields of the TLRs, which 
provides the feasibility of observing the strong gauge field 
effects, such as HB spectrum, in a realistic hybrid solid- 
state system. 

As illustrated in Fig. 1, the system we study is a 2D 



square lattice, where the basic unit consists of a NVE 
confined in a microwave-driven two-mode TLR with the 
length L, the inductance Ft and the capacitance Ct, a 
narrow center conductor and two nearby lateral ground 
planes. Each NVE has the average NV center density 
1 X 10^^ cm~^ ^34|, where each NV center is negatively 
charged with two unpaired electrons located at the va- 
cancy, which can be modeled as a two-level system in 
the ground-state subspaces as shown in Fig. 2. For 
clarity of description, we adopt following denotations: 
|^A,ms = O) = \g), and \^A,ms = +l) = |e) (we dis- 
tinguish the degenerate sublevels pA, ms = ±l) by ap- 
propriate external magnetic field or polarized irradia- 
tion). In our system, microwave-photons can tunnel be- 
tween adjacent TLRs via capacitive coupling by connect- 
ing TLRs via a central coupler, which serves as individual 
tunable quantum transducers to transfer photonic states 
between the resonators. 

The overall system is governed by the Hamiltonian 
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is for the {p,q)-th TLR's 
the creation (annihilation) oper- 



ators of the full-wave mode of the resonator and ujc — 
27r/(V-FtCf) the corresponding eigenfrequency. The col- 
lectively magnetic coupling between NVE and all the mi- 
crowave modes can be modeled as 
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where g'p ^ is the single NV center's vacuum Rabi fre- 
quency in the (p, q)-th. resonator for mode c, and is the 
frequency of the external driving laser with the Rabi fre- 
quency rip^g and the phase Op^q- The collective transition 

operators of NVE are defined as 5+^ = Yl!j=i \^)p,qp,q{9\ 
and Sp^q — J2f=i \9)p.qp,q{^\ '^ith Np^q the number of the 
NV centers in a NVE inside the {p, g)-th resonator. 



The Hamiltonian of a NVE reads H 



NV 



2^eJp^q 



with We the energy-level spacing of the states |e) and 
\g). All the spins in NVE interact symmetrically with a 
single mode of electromagnetic field by magnetic-dipole 
coupling because the mode wavelength is larger than the 
spatial dimension of the NVE if the spin ensemble is 
placed near the TLR's field antinode. We suppose that 
the mode Up^q couples with the vertical adjacent sites, 
and bp^q couples with the horizontal adjacent sites. So the 
tunneling between adjacent sites can be governed by the 



Hamiltonian H, 



(p,q) 



f^iTadp q^ldp^q + Tall 



Tiibp+i,qbj, q + Ti,bp^qbp_^_i q). Here the intercavity photonic 



tunneling rates Ta{Tt,) are tunable experimental parame- 
ters due to the flexibility of the central coupler. The cen- 
tral coupler may be conceived as a current-biased Joseph- 
son junction phase qubit [s^, or Josephson ring circuit 
[3], or a capacitive coupling element fspl . or an active 
non- reciprocal devices as proposed in |37| . 

Under the strong driving case flp^q ^ Tc, \/N^gp q, 
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we define the new energy levels as \+)p,q = (e*^"" \e)p^q + 
\g)p,q)/V2 and = (-e*^-«|e)p,g + \g)p,g)/V2. In- 

troducing Ac = Wc — We and Sp^^ = Ac — ^p,q, and sup- 
posing \Sp^q\ <C ^p,q, we have the effective Hamiltonian 
of the mode a in lattice site (p, q) under the new ba- 
sis as = Sp^^a^p^^ap q + {\gp^q& ^^'''''^p.q^p.q + H.C.), 
with S+q = Y!j=i l+)p,qp,g(-| (sGc SccI in supplemen- 
tary material). On the other hand, the effective Hamil- 
tonian describing the mode h is Hp^ = Sp ^bj^^bp g — 

For simplicity, we suppose 6p ^ — S, Tc — T, and 
9p,q — 9 for all p,q,c. In the limit of (5 ^ T, g, we 
can adiabatically eliminate the photonic modes and get 
the effective Hamiltonian between NVEs (see Sec. II in 
supplementary material), 

H,s = - J'{J2 e'^'--^-'-^'-^^S+Jp+i,q 

P,q (2) 

+ s+^Sp,q+i ^ jj 

where / = T{g/25f. 

Using the Holstein-Primakoff transformation 

= Bl q^ Np^q - BlqB p^q V^^.,' 

^p,q ^ Bp,q\J ^p,q ~ Bp,qBp,q — \J ^p,qBp^q, and 
S^^^ = {Bl^Bp^q - Np^q/2) [Hi, where the operator 
B^ ^{Bp^q) fulfills the bosonic commutation relation 
[Bp^q, Bp, ~ Spq^p'qi in the case of low number of NVE 
excitations, we can map the collective raising (lowering) 
operators 5'^^ of the {p,q)-i\i NVE into the bosonic 
operators ^ and Bp^q. These transformations change 
the Hamiltonian H^f / into a concise form as 

Hb = -J{Y,e'^'^-'~'^^'-'^BlgBp+i^q 
p,q 

+e''^'--+'-''--^Bl^Bp^^+^ + H.C.), (3) 

where J = Np^qJ = T{geff /2S)'^ and the collective cou- 
pling rates g^ff between NVE and TLR in each lattice 
site are supposed to be equal in our case. 

Because of the spatial variation of the tunneling phase, 
the wave function of a polariton from one lattice site to 
another acquires a nontrivial phase, which can be inter- 
preted as an effective Aharonov-Bohm phase. With a 
suitable tuning of the tunneling phases between pairs of 
neighboring NVEs, the polariton turns out to experience 
a non-trivial artificial gauge potential A, which can be 
identified by 

ep,q = - / AdL=- A-dL, (4) 

" JQ,q Jp,q 

where the integral is performed along the segment con- 
necting the TLRs, and e is the effective charge. With 
the choice of the symmetric gauge, we set Op^q = —irpqa, 
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FIG. 3: (Color online) For 10 x 10 lattice, the spatial distri- 
bution of absolute square of the polariton Bp^q wave function 
in the ground state of Hamiltonian (3) with open boundary 
conditions, under different values of a: (a) a = 0.181. (b) 
a = 0.182. 
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FIG. 4: (Color online) The refation between the critical value 
of ao and the system size A^. 



where a = $/$o with $ the magnetic flux through a 
unit cell, and $o = h/e the flux quantum. So a uni- 
form artificial gauge field will emerge during the polari- 
tonic dynamics. The sum of the tunneling phases along 
a closed loop surrounding the plaquette is 2{9p+i_q + 
dp,q+i — dp,q — 0p+i^q+i) = 2Tra, which is actually the 
flux quanta per plaquette, namely, the strength of the 
artificial gauge field. Note that the phase errors of the 
local microwave source only cause the local gauge field 
variance. The total fiux remains unchanged. 

Next we focus on the observable consequences of a uni- 
form magnetic field in the present system. For simplic- 
ity, we consider the case of non-interacting polaritons in 
finite-size 2D square lattice under the tight-binding limit 
and the single-polariton subspace. As shown in Fig. 3, 
we plot the spatial distribution of polaritons Bp q in the 
ground state of the 10 x 10 lattice under different val- 
ues of a. When a = 0, the Hamiltonian in Eq. (3) 
can be reduced to that for a free particle in an infinitely 
deep square potential well and the form of ground state 
takes the Sine functions. When a increases from to 
ao ~ 0.181, we find that the polaritons concentrate at 
the central region which implies the wave-function of the 
ground state gradually change from Sine function into 
that of the 0-th Landau level with zero angular momen- 
tum m, and the radius of the 0th Landau level decreases 
as the gauge field increases. This transition which is de- 
noted by a is really from 0-th Landau level of m = to 
0-th Landau level of to > 0. However the distribution 
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changes dramatically once the value of a exceeds a criti- 
cal value OfQ J where the 0-th Landau level with to > mo- 
mentum becomes ground state. The physical mechanism 
behind this intriguing phenomena is the competition be- 
tween the gauge field and the size eiTect. As shown in Fig. 
4, we find that the crossing point ao can be well fitted by 
the equation ao = 2/ (TV -I- 1) for the N x N systems with 
N > A (see Sec. Ill in supplementary material). If we fur- 
ther increase a, many level crossing points would appear. 
We find that all the crossing points are near certain ra- 
tional points a = pIq, at which the wave-functions show 
regular oscillations 2| . After Fourier transformation, one 
can see that such regular oscillations come from the co- 
herent interference between the peculiar points in mo- 
mentum space (see Sec. Ill in supplementary material). 
It is very interesting that the results from the model with 
the periodic condition in thermodynamic limit can be ob- 
served in this small and open boundary system. In order 
to measure the level crossing, we can couple the NVEs to 
superconducting qubits, and transfer the state of NVEs 
to the nearby qubits. Then we can measure the ground 
state population distribution of polaritons by measuring 
the states of superconducting qubits. 

To investigate the possibility of observing the fractal 
band structure of HB, we diagonalize the Hamiltonian 
Hb in the single-polariton subspace for several small lat- 
tice sizes. As shown in Fig. 5a, the HB structure can 
be observed in the size of 5 x 5. It becomes clearer in 
the whole transmission spectra with the growth of the 
size. Experimentally, transmission and reflection mea- 
surements have been by far performed routinely in small- 
scale circuit QED systems with one or two resonators 
(2ft . ,40.] . We can characterize our hybrid circuits by spec- 
troscopic measurements, and a feasible method is to mea- 
sure the resonator transmission through a network ana- 
lyzer. In realistic experiments, the transmission spec- 
trum of these superconducting resonators will show a se- 
ries of Lorentzian peaks, whose central frequencies cor- 
respond to the eigenfrequencies of the Hamiltonian Hf,f f 
in Eq. (3). Alternatively, using homodyne detection fol- 
lowed by sampling and averaging after amplification, we 
may reveal the coupled resonator-spin dynamics by mea- 
suring the amplitude of the exponentially damped mi- 
crowave signal that leaks out of the resonators after its 
interaction with the spin ensembles 4l| . So far it is still a 
great challenge to observe such a fascinating structure in 
ordinary solid-state system, where the major challenge 
comes from the extremely high magnetic field required 
in realistic experiments [42| . In this sense, our proposal 
opens an entirely new arena to investigate condensed 
matter and many-body system with light. Meanwhile, 
the experimental progress, especially in large-scale cir- 
cuit QED [3l|, HI] , has raised the possibility of observing 
the HB spectrum directly in realistic experiments [i^ . 

Finally, we survey the relevant experimental parame- 
ters. First, the TLR cavity with an inductance Ft = 45.6 
nH and capacitance Ct = 2 pF leads to a full wave fre- 
quency ujo/2tt = Dgs + Ac = 3.31 GHz with Ac/2tt = 440 




FIG. 5: (Color online) The energy spectrum for (a) 5 x 5 
lattice; (b) 6 x 6 lattice; (c) 8 x 8 lattice; (d) 10 x 10 lattice, 
where energy is in units of tunneling strength J. The vertical 
axis is the eigenenergy of the system, and the horizontal axis 
indicates magnetic flux through a unit cell. 



MHz. Second, our scheme requires the large-detuning 
condition, namely, the detuning Ac S> geff- Third, 
the classical field should be tuned to J7/27r — 400 MHz 
3> geff, which makes the energy shift induced by the 
cavity mode negligible f47]. So if we take the values 
of the parameters geff/2n = 8 MHz, T/27r ^ 4 MHz, 
and 6/2n = 40 MHz, yielding J/27r = Tgl^^/iS^ = 0.04 
MHz, our scheme is feasible. This could be confirmed 
by recent experimental demonstration of coherent cou- 



pling of a NVE and a superconducting flux qubit [48 1 
as well as the experimental advances in excellent quan- 
tum control with strong magnetic coupling (~ 27r x 10 
MHz) between TLR and NVE and the cavity linewidth 
of K/27r ^ kHz (22| - [25| . On the other hand, the electron 
spin relaxation time Ti of NV centers ranges from 6 ms 
at room temperature [i^ to 28 ~ 265 s at low tempera- 
ture [131 . In addition, the dephasing time T2 > 600 ^s for 
NVE with natural abundance of ^^C has been reported 
51 1 . A later experimental progress [52| with isotopically 
pure diamond sample has demonstrated a longer dephas- 
ing time to be T2 = 1.8 ms. Therefore, the parameter J 
is higher by nearly two orders of magnitude than the dis- 
sipation rates of the NVE, which makes reliable quantum 
simulation feasible. 

Compared with previous cavity QED protocols [Tsl - lssj . 
our scheme requires only one driving source at each site, 
and the effective couplings between neighboring spins are 
much larger than in Ref. because the excited states 
are not required to be adiabatically eliminated in our 
case. However, so far we have ignored the detrimental in- 
fluence from the nuclear spin, such as ^^C defects, in the 
NVE, which could be alleviated by isotopically purified 
^^C diamond through the purification technique |26l,[5^. 
Another decoherence source is the dipole-interaction be- 
tween the redundant Nitrogen spins and the NV centers, 
which could be reduced by improving the nitrogen to 
N-V conversion rate while maintaining the large coUec- 
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tive coupling constants [25|. Alternatively, this problem 
could be overcome by applying the external driving field 
to the electron spins on the Nitrogen atoms. It would 
increase the coherence time of the NVE if these spins are 
flipped on a time scale much faster than the flip-flop pro- 
cesses [2^. On the other hand, the dephasing time can 
be greatly enhanced by decoupling the electron spin from 
its local environment with a spin echo sequence. Using 
this technique, the dephasing time of the NVE reaches 
3.7 fj,s at room temperature 25]. 

In summary, we have discussed how to simulate gauge 
fleld in a J-C lattice of NVE-TLR system, where the 
motion of polaritons (collective excitations) in NVEs is 
analogous to the motion of charged bosonic particles in 
a magnetic field. We have discussed the competition be- 
tween the size effect and the gauge field on the ground 
state crossover, and also discussed the possibility of iden- 
tifying signatures of the HB in the optical spectra of 



the superconducting resonators. With currently avail- 
able techniques, we argue that our system lends itself as 
a well-suited quantum simulator for investigating phe- 
nomena encountered in condensed matter physics, and 
our study would be useful for the future spintronics tech- 
nology. 
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I. HAMILTONIAN FOR THE SINGLE-SITE NVE-TLR 



We suppose that each NVE couples to a two-mode TLR in each site. By applying an 
external resonant driving microwave on the NVE, the Hamiltonian of the (p, g)-th site reads 

c=a,b j=l j=l 

+ Zl^p,<?Cp.5l^)p-9P-'/(^l) + H-C.], (1) 

c=a,b 

where Qp^q is the Rabi frequency of the driving microwave, and 6p^q is the relative phase. 
In the rotating frame regarding Hq = T.c=a,b^ecl^q'^p,q + ^eYlj=i |e)^,g^,g(e|, we get the 
Hamiltonian as 

Hp,q = z2^^^~ ^^H,g^P'i + lz2^^^ \eyp,gp,q{9\ + 22 ^P,gCp,g|e)^,g^,9(^|) + H.C.]. 

c=a,b j=l c=a,b 

Here we define = uJc — uie for c = a, b. 

Under the strong driving case Qp^q ^ g^^, we can define the new energy levels as \+)p,q = 
(e*^f'''|e)p,g + \g)p,q)/\/2 and |-)p,g = (-e*^'''" |e)p,g + \g)p^q)/\/2. The Hamiltonian of the 
system under the new basis is 

Hp,q = J2 ^cCl^P. + + I - 

c=a,b j=l 

+1 E 9;,qCp,,e~^'''H\+yp,, - (-1) 

c=a,b 

+1 E + (-!)• (2) 

c=a,b 

Let's neglect the energy shift induced by the cavity mode, which is much smaller than the 
energy shift induced by the driving laser. Then we obtain 

Hp,q = ^cCpqCp^q + ^ ]{ - (|+)p,<jp,<j(+| ~ |~)p,<jp,g(~|) 

c=a,b j=l 

+ 2 E E [9p,qCp,q^-''''''i\+)lAq(-\ - + H.C.]. (3) 

j=l c=a,b 

Suppose Aa — \^p,q\, Ah ^ — |r2p_g|, and flp^q ^ \/ ^p,q9p,q^ obtain the effective 
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Hamiltonian under the rotating wave approximation, 

Zjeff = of o + ftt ^ ^ _\^\r,'' a ^ (-\ 

J = l 

-<A.e-^''"1+)U.(-| + H.C.], (4) 

where 5^^^ = - VLp^q and 5^^ = + fip^g. Let's define 5+^ = Xlj^i" and 
^p,q ~ Yl!j=i for (p, g)-th site. We denote the Hamiltonian of the couphng between 

mode a and NVE by iJ^g = Spgal^ap^g + (|fi'p,ge"*^P'''ap,g5'+g + H.C.). On the other hand, 
the effective Hamiltonian describing the coupling between mode b and NVE is given by 

-'-'p,g ^p,q^p,q^P,q \2^P,q ^p,q'^p,q ^ ) ■ 

In the rotating frame regarding Hq = ^p ^aj ^ap^g + Sp gbl gbp^g, we rewrite Eq.(4) as 

/T^// = ifl" f,-ii^p,i+^p.,qt)a ,9+ --o^ f,-iiSp,i-K,<i^)h ,9+ + HC (5) 

P>5 2^''^ "'P^I'^P^q 2 '^P,q'^p,q^ W/* 



II. EFFECTIVE HAMILTONIAN FOR COUPLED CAVITY SITES 

The coupling between [p, q)-th site and the nearest TLR sites can be modeled by following 
Hamiltonian. 

TV 

Ht = - ^ h{Taal g^^ap^g + Taa^p ^ap^g+i + Tkbp+i^gb^ g + Tbbp^gbl^^ g). (6) 

p,q=l 

We diagonalize the Hamiltonian (jH]) with the Fourier transformation a^^^p = 

^ Z^q=l ^ "P>9 '^^^'-^ "kp,q ^ Z^p=l ^ '^P,qy 

2tv N 2n N 

Ht = - Y.'^'^A^qKi - Yl 5Z^fcp4„p«fcp,p> (7) 

kp=2n/N q=l kp=2n/N p=l 

where T^^ = 2r;, cos(27rA;p/A^) and Tk^ = 2TaCos{27ikg/N). We define Fp^kp = 6-''^^''^^'^ and 

P _ p-i2TTkqq/N 
J^q,kq ~ 

It is found that the coupling between vertical sites and between horizonal sites can be 
calculated independently. The effective Hamiltonian of the system for the vertical line q is 

N 2-K 

H. = - E(2^J-^^"^^'^'^"'^"*^^^'^^^"^ + H-^-) - E Tkpbl,Kr (8) 

p=l kp=2TT/N 
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For simplicity, we assume 5^ ^ = ^ = S, g^ ^ = g^ ^ = g, and Tp^g = T. With the transfor- 
mation 6fcp,<; = Y.p=i Fk^,pbp,q, we get 

N 2tx ^ 
p=l kp=2n/N ^ V iV 

In the hmit of 5 ^ T, g, we can adiabatically ehminate the cavity modes b^p^q and get 

p,p' kp=2n/N 

In the limit of 5 > T, we get 1/(5 + T^J ^ ^ - ^, Efc, ^fc^.P^^y = and 

Y,kpFkp,pF^^^p/Tkp = T((5p,p'+i + 5p+i,p'). Neglecting the Stark shift terms, we get the ef- 
fective Hamiltonian as 

^v" = E -^^^'^'''^'''''''''^^UiJp,, + H-C- (10) 



Similarly, we can calculate the effective Hamiltonian for the horizontal line p by adiabat- 
ically eliminating the cavity modes as the following form. 



Combining Eqs. ([8]) and f lTT]) . we obtain the effective Hamiltonian for the whole system as 

i^eff = - /(Ee^^''"'"'''"''^'^pt.^;+i,, 

P.9 (12) 

+ e^^^-+-^-)5+,5,^,_,i + H.C.), 



where j' = T(g/25) 
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III. THE GROUND STATE NEAR LEVEL CROSSING 

We find that the ground state level crossing can be observed in the system as small as 
5x5 lattices. As shown in FiglH the critical value of gauge field is ao = 0.333 ^ 1/3. 
The numerical result of the ground state near the critical gauge field ao for the 5x5 
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FIG. 1: (Color online) The spatial distribution of the NVE polaritons for a 5 x 5 lattice 
under different values of a. (a) a = 0.334. (b) a = 0.333. 

lattice system is shown as follows, 

$(0,333) =0|l)i,i + (0.0003 - 0.0002z)|l)i,2 + (0.144 - 0.0009z)|l)i,3 
+ (0.0003 + 0.0002z)|l)i,4 + 0|l)i,5 

+ (0.0003 + 0.0002i)|l)2,i + 0.1448|1)2,2 + (0.3536 - 0.0011i)|l)2,3 

+ (0.1448 - 0.0009i)|l)2,4 + (0.0003 - 0.0002i)|l)2,5 

+ (0.1440 + 0.0009z)|l)3,i + (0.3536 + 0.0011z)|l)3,2 + 0.5772|1)3,3 

+ (0.3536 - 0.001H)|1)3,4 + (0.144 + O.OOOQz) 11)3,5 

+ (0.0003 - 0.0002z)|l)4,i + (0.1448 + 0.0009z)|l)4,2 

+ (0.3536 + 0.0011i)|l)4,3 + 0.1448|1)4,4 + (0.0003 + 0.0002i)|l)4,5 

+ 0|l)5,i + (0.0003 + 0.00022)11)5,2 + (0.1440 + 0.0009i)|l)5,3 

+ (0.0003 - 0.0002i)|l)5,4 + 0|1)5,5. 
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$(0,334) = - 0.0004|l)i,i + (0.1251 - 0.00713i)|l)i,2 + (0.0022 - 0.1763i)|l)i,3 
+ (-0.1232 + 0.00745z)|l)i,4 + 0.0004|l)i,5 

+ (0.1251 + 0.0713z)|l)2,i + 0.3065|1)2,2 + (0.0018 - 0.2890i)|l)2,3 

+ (0.3064 - 0.00392)11)2,4 + (0.1260 - 0.0698z)|l)2,5 

+ (0.0022 + 0.1763i)|l)3 1 + (0.0018 + 0.2890i)|l)3 2 + 0|1)3 3 

(14) 

+ (-0.0018 + 0.2890i)|l)3,4 + (-0.0022 + 0.1763i)|l)3,5 
+ (-0.1232 + 0.0745i)|l)4,i + (-0.3064 + 0.0039i)|l)4,2 
+ (-0.0018 - 0.2890z)|l)4,3 + 0.3065|1)4,4 + (0.1242 + 0.00729z) 11)4,5 
+ 0.0004|l)5,i + (-0.1260 + 0.0698z)|l)5,2 + (-0.0022 + 0.1763z)|l)5,3 
+ (0.1242 - 0.0729i)|l)5,4 - 0.0004|1)5,5. 
Here |\I/(a)) denotes the ground state at the field a, and |l)p,g denotes the Fock state |1) 
of the polariton at the lattice site p, g . It is easy to verify that (<l'(0.333)|$(0.334)) = 
1.58 X 10~^ ~ 0. Therefore, there is level crossing near the critical gauge field = 0.333. 

In general, we can define the ground state fidelity for the system F{a) = (<I>(a)|$(a + 5)), 
with 5 — 7- 0. For those F{a) approach zero, the level crossing occurs. We have numerically 
calculated the F{a) and identified the level crossing points for different sizes of lattices. For 
example, for 5 x 5 lattice, the level crossing point is near 1/3. For 6x6 lattice, the level 
crossing points are near 2/7, 3/8, 2/5 and 1/2. The critical points are fitted very well with 
rational number p/q for different lattices, and the first crossing point is just 2/(L + 1) 
for the L X L systems (L > 4). To qualitatively learn the nature of the first level crossing 
at ao = 2/(L + 1), we estimate by perturbation approach. 

First we consider a zero magnetic field case, a = 0. Now the effective model in Eq.(12) 
is reduced to a continuum form as 

where n^' = jj- Considering the boundary condition, 

^(0,0) = 0, ^(Lp,0) = 0, 
^'(0,L,) = 0, ^(Lp,L,) = 0, 
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where Lp (Lg) is the number sites along p(q) direction. We have the wave-function as 



2 . fT^p^ . fT^q 

and the energy eigenvalues as 



^(p,<?) = ^==sin(— )sin(— ), 



where p = 1,2,3, ..Lp and q = 1,2,3, ..Lg. It means that the ground state for the a — >■ 
case is labeled by p = 1, g = 1 with the energy 

The first excited state is labeled byp = 2, g=lorp=l,g = 2 with the energy 
mm(i?(,,),i?(,,))=mm(-[(— ) +(77)])- 

Secondly we add a tiny magnetic field a. The energy spectrum of the system changes 
little. So we may use the perturbation theory to deal with this case. Here we choose the 
Landau gauge, A = {By, 0), where the gauge field strength B = 2Txa. In the following part 
we set = Lx = L. 

Now the model in Eq.(12) in a continuum form turns into 

H'r. ^ — + (16) 

^ 2m 2m ^ ^ 

In the perturbation theory, we use the plane waves of non-magnetic case ^I/(p,g) to describe 
above model. Then for the lowest eigenstate ^(1,1), and the excited state "^{2,1), we have the 
corresponding perturbation energies as 



9=1 



term, 



and 

1 r / . On 1 r / 27r , On 1 1 X ^ 2 TT , , O 



9=1 

„2 



where the a terms can be neglected when a ^ 1. So there may exist level crossing between 
and \I/(2,i) when 

-^(1,1) = -^(2,1) 



a=l/2 




that is 



FIG. 2: The total fermion density of vr-flux case on 20 x 20 lattice. 



1 r,27r,9, 1 r.n.o^ 1 1 sr-^ 27r , , 1 1 , 

9=1 <?=! 



for tiny a. That is 



3 TT 2 1 1 



2m^L' L2 2m5^'Z^'™^^) 

9=1 



or 



3 1 



ar 



2(L + 1) 



with X]g=i ? = Therefore we estimate the first level crossing at = 1.5/(L + 1) by 

perturbation approach, which is qualitatively same as the numerical result = 2/{L + 1). 



IV. INTERFERENCE PATTERN FOR a = p/q 

From the ground state fidelity we have found that all the level crossing points are near 
certain rational points a = p/q, at which the wave-functions show regular oscillations. 
Fig. 2 illustrates the particle density of a = | case on 20 x 20 lattice. After Fourier 
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FIG. 3: The illustration of interference pattern of vr-flux case on 20 x 20 lattice in 
momentum space. The blue lines are the boundary of reduced BZ. The 8 resonant points is 



transformation, we find that such regular oscillations come from the coherent interference 



on an L X L lattice with open boundary condition, the particle density can be transformed 
into 8 resonant points inside the reduced Brillouin zone (BZ) of momentum space as 



around the nodal points {kp, kg) = (|, ±f ). Let's explain this result. For an infinite system. 



the L X L lattice with open boundary condition as an infinite system with periodic L x L 
super-cells, of which the wave-functions are fixed to be zero at the interfaces between two 
Lx L super-cells. Thus, modulated hj Lx L super-cells, each nodal points split to 4 points. 
As a result, we can read out the information of a model with the periodic condition in 
thermodynamic limit in this small and open boundary system. As shown in Fig. 3, for a 
20 X 20 lattice, the 8 resonant points in the reduced BZ zone are already very close to the 
2 nodal points for an infinite system. 




between the peculiar points in momentum space. Thus for the ground state of a = ^ case 




the case of a = | is just a vr-flux model, of which there are two nodal points in the reduced 
BZ as (|, ±1 ). To transform the wave-function of the system with finite size, we can consider 
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